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Abstract 

Relativistic plasma with radiation at thermodynamic equilibrium is a general sys- 
tem of interest in astrophysics and high energy physics. We develop a new self- 
consistent quasi-particle model for such a system to take account of collective behaviour 
of plasma and thermodynamic properties are derived. It is applied to electrodynamic 
plasma and quark gluon plasma and compared with existing results. 



PACS Nos : 05.30.-d, 52.25.Kn, 12.38.Mh, 52.40.Db, 12.38.Gc, 05.70. Ce 

Keywords : Relativistic plasma, thermodynamics, quark gluon plasma, quasi-particle. 







1 Introduction 



A typical clcctrodynamic plasma of interest in astrophysics [1], consists of electron, positron 
and photons (e, e"*", 7) at thermodynamic equilibrium. Since it is a plasma medium, individ- 
ual particle properties are modified by collective effects of plasma. One way to take account 
of this effect, is to use quasi-particle model. Just like in Debye theory of specific heat or 
theory of liquid helium etc. thermal properties of the medium may be viewed as a result 
of thermal excitations or quasi-particles, like plasmons and dressed photons as a result of 
quantization of plasma waves and electromagnetic waves in plasma. The standard procedure 
is to obtain the classical dispersion relations for plasma or electromagnetic waves in plasma 
and on their quantization, we get quasi-particles, namely quasi-fermions corresponding to e, 
e+ and quasi-bosons like dressed photons. We study the statistical mechanics and thermo- 
dynamics of such a system of quasi-particles. One such study was attempted by Medvedev 
[1] which we modify and correct it to get the present model and then generalize it to quark 
gluon plasma (QGP). 

QGP is a plasma made up of quarks and gluons [2], governed by strong interaction 
called quantum chromodynamics (QCD). It is similar to (e, e"^, 7) system with electrons 
(positrons) are replaced by quarks (antiquarks) and photons by gluons. QED (quantum 
electrodynamics) is replaced by QCD. QGP also exhibits collective behaviour and hence the 
thermodynamic properties of QGP are modified and may be studied using quasi-particle 
model. Of course, this kind of study was first attempted by Peshier et. al. [3] and latter 
it was modified with corrections and studied various groups [4]. However, initially, these 
models were thermodynamically incosistent [5] and also not self-consistent in calculating the 
thermodynamics (TD) and quasi-particle properties. Few attempts were done to correct 
these inconsistency [5, 4] by reformulating statistical mechanics (SM). However, in Ref. [6], 
we have developed TD consistent model of QGP without reformulating SM. Here, in this 
new model, all previous problems of quasi-particle model, namely TD inconsistency and 
self-consistency, are solved. 
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2 Quasi-particle model of plasma 



Both the systems we discussed are highly relativistic and hence we develop here quasi- 
particle model for such ultra-relativistic systems. We assume that the collective excitations 
of plasma leads to a system of non-interacting quasi-fermions and quasi-bosons, obeying 
Fermi and Bose statistics respectively. Following the standard statistical mechanics [7], 
density of quasi-particles may be written as, 

where z is the fugacity, =F refers to Bosons and Fermions. is the degenarcy associated 
with the internal degrees of freedom. Here is the energy of quasi-particle which may be 
obtained from classical approximate dispersion relation as, 

efe = hsJ¥J+Jl , (2) 

for Bosons and 



ek = h,J¥J+^ , (3) 

for Fermions. These forms of energy-momentum relations are widely used in quasi-particle 
models of QGP [3, 4] with uj^ replaced by temperature dependent masses, which they obtain 
from the finite temperature field theory calculations in ideal thermal bath. The general 
expressions for are very complicated even at high momentum limit and following Medvedev 
[1], we approximate them to above simpler equations with an error of about 3%. c is the 
speed of light and h is the Planck constant, ujp is the plasma frequency, given by. 



= ^ = a j7 , (4) 



for (e, e"*", 7) system, e is the charge and rig the electron density which is also equal to 
positron density for a system with chemical potential equal to zero. In quasi-particle models, 
by definition, electron density is same as that of quasi-electrons. T is the temperature of the 
system. 
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3 (e, e+, 7) system: 



Let us first consider (e, e"*", 7) system with chemical potential zero, or 2; = 1, and the density 
of quasi-electrons is, 

n 1 



o / t^'^'^^ , , (5) 

which may be rewritten as 

ne^l^i^r Hdxx^ , ^ . (6) 



e V +1 

The density of positrons is same as that of electrons. Similarly, the density of quasi-photons 
may be written as, 

n^^p-{^f rdxx'' , ^ , (7) 

^ 27r2 ^Uc' Jo J^2^an-'^ _^ ^ ' 

where the fugacity is one for photons. These equations need to be solved self-consistently 
because Ue which is to be determined is inside the integral through ujp. Redefining the 

variables, the final equation to be solved self-consistently is. 



/■oo 1 00 f — iV'-l) 

fe-l dxx' ^a^flY.^-^K,{alU), (8) 



where 



-2 _ 4^e 

a = — — a 



and 

2 _ 2 TT^ {h c)^ Ue 

a is the fine structure constant and K is the modified Bessel function. Once we know Ue or 
/g , we can obtain photon density from the relation, 

n,-^Af/d^fet]Malfe), (9) 
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which follows from the Eq. (7). Similarly, energy densities are given by, 

rp4 



27r2 {hcf 



or 



rp4 



Jo 



dxx 



ifjfK,{afel) + 3 {afelfK2{afel)]) 



in terms of modified Bessel functions K, for electrons and 

rp4 



(10) 



(11) 



or 



^ 

27r2 {hcf 



Jo 



dxx 



x^ + /g 



OO -1 

^- [(a/eO'i^i(a/eO + 3(a/e/)'i^2(a/eO]) , 



(12) 



(13) 



^ 2 7r2 (^c)3;e^ 

for photons. The blackbody Planck's distribution may be easily read from above equations 
as 

rp4 



de^{x) — 



9i 



27r2 (he) 



X 



(14) 



where x = huj/T. Thus the Planck's distribution is modified for the non-zero value of a due 
to plasma with a cutoff frequency related to plasma frequency. 
On taking a limit a — > in Eq. (6, 7), we get 



and 



where C(3) is a Riemann zeta function and has a value approximately 1.2 and r){3) is related 
to C(3) as r](3) = | C(3)- Hence the Ue and are related expected for ideal 

gas. Note that in the earlier calculations, [1], the starting point of the formalism is ng = | n^, 
which is not right and such a relation is for energy densities. Secondly, it is an ideal gas 
relation which need not be true for non-ideal system that we are discussing. Here in our 



model we don't use such relations between Ue and n^, rather both of them are calculated 
self-consistently by solving the coupled integral equations, Eq. (6, 7), 

Other thermodynamic functions, like pressure may be obtained from the thermodynamic 
relation 

dP 

e-T—-P, (15) 

which on integration gives P = ^e, same as that of ideal relativistic gas. Again this result 
also differs from Ref. [1] where their expression for pressure is not valid for massive particles. 
It is interesting to look at the expression for plasma frequency, 

u^l-:^9eafAlr^^-^e^T\ (16) 

where we took /g ^ 2 77(8) = |C(3). It is very close to |e^T^, obtained from the finite 
temperature field theory calculations, but note that it is not exactly equal because of our 
self-consistent calculations. 

4 Quark gluon plasma: 

Prom above discussion on (e, e"*", 7) system, it is clear, in contrary to earlier work [3, 4, 5], 

one need to study QGP self-consistently in quasi-particle model. QGP is expected to form 

in ultra-relativistic heavy ion collisions and may exist in astropysical objects like neutron 

star, early universe etc. [2]. The lattice gauge theory (LGT) simulation of QCD predicts 

such matter, but seems to be non-ideal plasma near to the transition temperature from 

hadrons to QGP. To explain this behaviour many phenomenological models, based on QCD 

and plasma, have been proposed and quasi-particle model is one of them. As we discussed in 

the introduction, many of quasi-particle models [3, 4, 5] lack self-consistency in calculating 

the quasi-particle mass and the thermodynamic properties. 

Let us consider, as an example, gluon plasma. The expression for gluon density, which 

is to be solved self-consistently, reduces to 
00 -1 

l = a'^-K,{alfg) , (17) 
1=1 ' 



where fg is similar to discussed earlier. Once we know fg, energy density may be evaluated 
from the relation, 

^p=^E74 [(a/./m(a/p0 + 3(a/,/)^i^2(a/,0]) , (18) 

expressed in natural units and gg — 16 for gluons. Note that, in QGP, gluons also carry 
charge and hence exhibits plasma oscillations even in the absence of quarks. Therefore, in 
pure gluon plasma, plasma frequency depends on gluon density. In QGP, it depends both on 
quark density and gluon density. Hence for gluon plasma we take Up — a Ug/T with a = a^ag 
such that a? — . otg is the QCD running coupling constant. Oq is a constant, equal 

to 87r/3 for (e, e"*", 7) system, but here we determine it by demanding that a;^ — > | g^T^ 
as T ^ 00, a pcrturbative finite temperature field theory result. So we get oq ~ 2.15. For 
constant a or a^, Sg oc and hence an ideal equation of state (EoS). However, LGT results 
show a modification of law near the transition temperature Tc and they actually used 
temperature dependent 2-loop order q;s(T'). Hence, we also consider similar expression to 
model q;s(T'), given by, 

<rj.. _ 67r / 3(153 - 19n^) ln(2 ln(T/A^)) \ 

^ (33 - 2n/) ln(T/AT) V (33 - 2n/)2 MT/At) )' ^ ^ 

where At is a parameter related to QCD scale parameter, nj, the number of flavors, is zero 
in pure gluon plasma to be compared with LGT results. At is the only parameter of our 
model to be found by fitting with LGT result. Once we know Sg{T), it is stright forward to 
get pressure P{T) by integrating the thermodynamic relation Eq. (15) and then all other 
thermodynamic quantities may be obtained. 

5 Results : 

For (e, e'^, 7) system we recalculated various thermodynamical quantities, reported in Ref. 
[1], using our model with proper corrections. The departure of these quantities from that of 
ideal system is too small to be noticed and tabulated in Table 1. Planck's distribution in 
plasma is plotted in Fig. 1 along with the Planck's distribution in the absence of plasma. 
Qualitatively similar results as in Ref. [1] with cutoff in the distribution etc. But, note that, 
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the numerical values differ and our values listed in Table 1 with plasma are smaller than 
without plasma. Where as in Ref. [1] just the opposite, seems to be a numerical error. 

When we apply our self-consistent quasi-particle model to QGP, as an exmple gluon 
plasma, we can explain very nicely lattice gauge theory (LGT) [8] results as shown in Fig. 2, 
using a single parameter to = ^t/Tc = -82. The self-consistently calculated plasma frequency 
is T"^ aX T = 5Tc, instead of ^g'^T^ used in all other quasi-particle models [4]. It is 

interesting to notice that it is close to the value we used in our earlier thermodynamically 
consistent quasi-particle model [6], \g'^T'^. 

6 Conclusions : 

We formulated a new self-consistent quasi-particle model to describe the thermodynamics 
(TD) of relativistic plasma, hke (e, e+, 7) and QGP. Basic idea is that because of the 
collective behaviour of plasma, TD of such a system may obtained by studying the TD of 
quasi-particles which are thermally excited quanta of plasma and electromagnetic waves in 
plasma. This is equivalent to a system of bosons and fermions with mass propotional to 
plasma frequency. Plasma frequency depends on the density, a TD quantity, which we want 
to find out and hence a self-consistent problem to be solved. Note that, in all earlier quasi- 
particle models of QGP [3, 4], this self-consistency is not taken into account. In addition 
they have problems of TD inconsistencies etc. and as a result of both problems, they need 
more than two parameters to fit the LGT results. Further extension of this model to QGP 
with finite number of quark fiavors, with and without chemical potential, may be future 
interesting problem. 

As long as (e, e"^, 7) system, studied earlier by Medvedev [1], is concerned, we modified 
his result with corrections in our new approach and obtained qualitatively the same results, 
but the numerical values and the sign of the modifications due to plasma are different. 
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Figure 1: The normalized Planck's distribution u{x) with plasma (continuous line) and 
without plasma (dashed line) as a function oi x = fiuj /T . 

Figure 2: Plots of s/T^ as a function of T/T^ from our model and lattice results (symbols) 
for gluon plasma. 



Table 1: Various thermodynamic quantities of (e, e"^, 7) system from our model with plasma 
and without plasma. 



Table 1 



Medium 




Ue [{T/hcf] 


n, [{T/hcf] 


Se [T7(M'] 


e, [TV{hcf] 


Plasma 


1.79925 


0.18230 


0.24164 


0.57527 


0.65705 


Vacuum 


1.80309 


0.18254 


0.24339 


0.57526 


0.65744 
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